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Abstract 



O , We consider flux lines in presence of a dense square array of columnar 

defects, with insulating core and radius of the order of the coherence length. 
The properties of the phase transition from the flux line regime to the local- 



> 

o. 

^^ , ized superconducting state are determined here. We show that the localized 

o 



superconducting state is a new vortex state because of the counterclockwise 
and of the clockwise supercurrent flows found around each defect. We propose 
here a general operator to count the vorticity of any cavity inside a super- 



'^ I conductor and give a new method to obtain the reversible magnetization of 

Q ', superconductors with multiply connected geometries. 



PACS numbers: 74.20.De 74.25.Dw 74.60. Ge 



Recent improvements in the resolution of electron beam lithography have allowed the 
fabrication of ordered arrays of columnar defects (CD) with radius R not much larger than 
the coherence length .^IteI. Flux lines (FL) in presence of a regular lattice of CD show new 
collective properties such as matching and multiple quantum trapping. Matching occurs 
every time the artificial CD and FL lattices are commensurate. According to some theoriegj, 
CD can trap many FLQ up to the saturation limit of R/2^(T), where T is the temperature. 

A regular array of CD turns the superconductor into a multiply connected geometry. This 
system presents the interesting property of surface superconductivity around the CD edges, 
which is crucially dependent on the properties of the medium exterior to the superconductor, 
as found long ago by Saint- James and DeGennea3. This medium, played here by the CD core, 
must be insulating such that the supercurrent flow is constrained to the superconductor's 
interioiill. For a non-perforated sample above the upper critical field Hc2 superconductivity 
disappears in the bulk. However in case of a perforated sample, such as the case of a regular 
array of CD, superconductivity remains above Hc2 around the CD edges, in the form of a 
thin layer of thickness ^(T), and this lasts until the field H^s is reached. Bezryadin, Buzdin 
and PannetierEl(BBP) determined several properties of this localized superconducting state 
(LSS), such as the ratio Hcz/ Hf.2 whose large R limit recovers the single planar interface 
problem, first considered by Saint- James and DeGennes. BBP also found that there are 
many possible LSS, each characterized by P, the number of FL trapped inside the CD. They 
determined that the transition between these LSS yields discontinuities in the magnetization 
(first order transition). Since they assume the LSS around each CD from the beginning, 
their conclusions only apply for fields above the transition. The interaction between the FL 
lattice and a single CD was also investigated below the transition where the local symmetry 
group caused by the deformation was determined^. 

In this letter we obtain the magnetic properties of the transition to the LSS, and propose 
a Ginzburg-Landau free energy expansion that describes the interaction between FL and 
a dense regular lattice of CD parallel to them with radius size typically of the order of 
the coherence length. The method also describes FL trapping by the CD, which shows 



a saturation inaxiinuni. In particular we determine the magnetic properties (Fig||T],0) of 
the FL matter in both sides of the transition, which is signaled in the magnetization by a 
slope discontinuity (second order transition). We propose a quantity, hereafter called fmh 
(Eq.(j^)), to count the number of FL trapped inside a CD for any B and T. Fig[0 shows that 
finh also determines the saturation number of FL inside the CD. finh is a spatial average 
value over the order parameter and of its gradient, that picks its non-vanishing contribution 
around the CD edge. In Figf^ we show the strip where the LSS is found in the diagram of 
the magnetic induction B vs. T. 

All results are obtained from numerical solutions of the present Ginzburg-Landau free 
energy expansion using the Simulated Annealing methodtd, which was first applied to 
superconductivity some time agdliJ. The magnetic properties were calculated using the 
Virial TheoremllHllla, whose original formulation already included the treatment of spatial 
inhomogeneitiealHl. The present treatment is for a regular square lattice of CD of side L, 
with a number A^ of FL inside the cell. FL are straight, thus there is translational invariance 
along the z-axis, and the problem is truly two-dimensional. The magnetic induction becomes 
B = {l/V)Jdv V X A = z N^o/L'^- We chose to work with a single CD {R = 2.0^o) per 
unit cell (L = 20^o); located at its center. This choice of unit cell does not allow the study 
of matching fields, because there is always commensurability between the two lattice param- 
eters. We work in the approximation of no screening, rendering this study interesting for 
an extremely high k, conventional superconductor such that London's penetration depth is 
much larger than the distance between two consecutive CD. 

The onset of the LSS is shown in FigQ], which displays the local supercharge and the 
local kinetic energy densities. Here B is fixed (A^ = 24), and two temperature are considered, 
T/Tc = 0.615 and 0.630, below and very near the transition, respectively. The transition 
occurs at T/T^. = 0.631. For both temperatures {ipl'^ is highly concentrated in the vicinity 
C,(T) around the CD edge, showing the onset of the LSS. For T/Tc = 0.615 small depressions 
of IV'P can still be identified outside the localized state around the CD. These depressions 
do correspond to FL external to the CD, because they coincide with local phase singularities 



and have counterclockwise supercurrent flow around them; these last two properties are the 
true signatures for the presence of a FL. In fact the high symmetry on the arrangement of 
external FL, formed around the CD below the transition, is a demand of the Saint- James 
and DeGennes condition. Clearly, symmetry of the external FL configuration around the 
CD implies that the total supercurrent vanishes at the CD centeiS. The supercurrent flow 
is best understood through the local kinetic energy density. For T /Tc = 0.615 there is an 
internal ring and several spikes around the CD, forming a highly symmetrical arrangement. 
For T/Tc = 0.630 the internal ring remains, but the kinetic density spikes collapse into a 
second ring. We have determined that this external kinetic ring is caused by a clockwise 
supercurrent flow, thus opposite to the flow of the inner original ring, whose counterclockwise 
flow is due to the P trapped FL inside the CD. Notice that the two kinetic energy density 
rings are located at the inner and the outer sides of I'i/'p in the LSS. 

We describe the superconductor with a square lattice of CD through the free energy 
density expansion f = F - B^/Sn = (l/V) J dv{ h^\Dtp\^/2m{x) + ao[T - T^{x)]\ij\^ + 
P\ip\y2 }, D = W - {2ni/^o)A, T^{x) = [1 - t{x)]T^, and l/m(f) = [1 - T{x)]/m. The 
CD are described in this model through the r(x) function, x being the distance from a CD 
center: r(x) = 1 inside and r(x) = outside the defect, respectively. To implement this 
numerical study, the following function is introduced, 

r(x)=2/(l+exp[(|x|/i?f]), (1) 

where i^ is a phenomenological parameter, such that the Saint- James and DeGennesQ con- 
dition corresponds to the limit K -^ oo. For the present numerical study we take K = 'b. 
It is convenient to introduce reduced units: the order parameter becomes A = 



/3/q;o(Tc — T)-?/;; both magnetic induction [B') and magnetization (M') are in units 
of ^/2HciT), Hc(T) /Atc = al{Tc — T)"^ / (3. The average energy densities are in units 
of Hc^{T)/At:: f, /^j^, and /^„^. To eliminate unnecessary temperature dependence, 
the following dimensionless fields are introduced: B = {1 — T /T^B' = (27c N/L'^)^qK, 
K = %/2TxV2^{TfH,{T), M = (1 - T/T,)M\ /,„, = (1 - T/T,)'fi^, etc. In reduced 



units the free energy becomes, 

/' = fL + ^Jdv { -[1 -r(f)/(l -T/TJ ]|Ap + |A|72 } (2) 

fL = l^Jdv[l-r{x)]aTf\DA\', (3) 



where ^(T) = ^o/Jl — T/Tc, ^o = yh /2maoTc. The corresponding Euler-Lagrange equa- 
tion is 

^{Tf (-Vr) ■ DA = -(1 - T)^{TfD^A - [1 - r/(l - T/Te)]A + |ApA. (4) 

For K -^ oo, r approaches the step function and so, in this hniit, Vr -^ —{R/\x\)6{\x\ —R)f. 
A physical solution of this equation only exists if its left side is finite and this demands that 
f ■ DA\\g\=fi = 0, which is the Saint- James and DeGennes boundary condition. This shows 
that the present theory provides a method to implement this boundary condition at the free 
energy level (Eq.(|])). The Virial Theorem, known to provide a simple method to understand 
Abrikosov's theoryo, gives the magnetization in this approximation of no screening, 

AnM' = i-fU^ + fiJA)/B' (5) 

fLH = ^ J dvi-x-VrMTflDAl' -\A\y{l-T/T^)], (6) 

where A is a solution of Eq.(^). //„^, the counter of trapped FL, directly contributes to 
the magnetization, thus explaining why the transition is of first order when the number of 
trapped FL changes inside the CD. The step function limit clearly shows that the major 
contribution to //„;, (Eq.(H)) comes just from the supercharge and gradient densities close 
to the CD edge. 

To find the solution of Eq.(P, first Eq.(^ is discretized, and then its minimum nu- 
merically sought through the Simulated Annealing methodlill. There are J values of the 
order parameter to be determined: A{nx, riy), n^ = 1, . . . , J, and n^ = 1, . . . , J, where the 
new working parameter, a = L/J^, is the distance between two consecutive points. All the 
present numerical results were obtained using a J = 30 array of points, so that a = (2/3).^o"- 

Our major results are displayed in the Figures. Fig|I| shows M vs. B for T = 0, 
corresponding to A^ = 1, . . . , 85 FL in the unit cell, and the inset shows the transition 



region. At the transition to the LSS (A^ = 68) an abrupt change of slope in the magnetization 
takes place. The discontinuities in the magnetization, signaled by Bp, P = 2,3,4,5, are 
the transitions to new trapped states, characterized by P FL inside. The B^ transition 
falls above the LSS transition, showing that multiple FL trapping and the LSS are two 
independent transitions. In fact the B^ transition has been utterly discussed by BBPS. The 
Bp transitions are best understood through Fig[0], which shows the diagram fmh vs. B. 
Notice that this curve naturally splits into five distinct lines, each associated to a P trapped 
state. For instance, for A^ = 1, . . . , 5 FL in the unit cell, and also for A^ = 9 one has that 
P = 1. But for A^ = 6, 7, 8, there are P = 2 trapped FL inside the CD. We believe that 
re-entrant effects are caused by the numerical procedure, which cannot resolve the very fine 
splitting between two possible states very close in energy, but distinct in the number of 
trapped FL. This difficulty is well resolved for B^ and -B5, but not for B2 and B4. The inset 
of Figpl shows finh vs. B for a higher temperature (T/Tc = 0.6), and is quite similar to 
the T = plot, except for the fact that the destruction of superconductivity happens at a 
lower field {B/k, ^ 0.5), which limits the highest trapped state to P = 3. Fig|0] shows that 
the LSS is found between the two data sets, which were obtained by keeping T fixed and 
varying A^. Above T/T^ = 0.8 there is no more LSS since CD interact collectively. The inset 
of FigP] shows the same discontinuity on the magnetization slope as the inset of Fig[|l|, but 
now obtained by varying T and keeping A^ fixed instead. 

Very near to a CD edge, the supercurrent fiow is counterclockwise, both below the transi- 
tion and also in the LSS. Remarkably, in the LSS, away from the CD edge, the supercurrent 
fiow fiips direction and becomes clockwise. Thus it vanishes at a critical radius |xo| from the 
CD center. Our numerical data is in agreement with a simple theoretical argument, given 
below, that this critical radius is 



\xo\ = ,JP/7rNL. (7) 

Below the transition consider the distance from the CD center to the symmetrical FL ar- 
rangement around the CD. External FL are attracted to the CD center, because CD are 

6 



pinning centers, and also repelled from them, because of the multi FL already trapped there, 
resulting on an equilibrium distance. This equilibrium occasionally can change when a new 
FL enters the CD interior. This process does not often occur because of the saturation 
limit, which is low for a CD with R ~ ,^o- Once the absorption process is completed, a new 
equilibrium symmetrical configuration of external FL is established. We find that when the 
transition line to the LSS is approached this equilibrium distance grows, meaning that near 
the CD edge the repulsion caused by the trapped FL becomes stronger and the pinning 
attraction to the CD weaker. The collective clockwise supercurrent flow arises between this 
symmetrical FL configuration and the CD edge. Within the LSS, and under the assumption 
that all localized states are independent, the wave function around each CD is approximately 
given hj ip = \ilj{\x\)\exp{iP6), where \x\ is the distance to the CD center, and P FL are 
trapped inside. Under the assumption of no screening in the unit cell {A^j. = 0,Ay = Bx), the 
supercurrent associate to this state, J oc |'?/'(|x|)p(P — it\x\'^B/^o)9, flips sign at a distance 
given by Eq.(|^. Since |a;o| coincides with the peak of \ip\'^, it must be of the order of C,{T) 
away from the CD edge. 

In summary we have studied here a dense FL lattice in presence of a dense square lattice 
of CD in the limit of no screening. We have obtained the magnetic properties of the transition 
to surface superconductivity around the CD. The interaction between the saturated CD and 
the external FL, the eventual capture of a FL by the CD are all described here through a 
Ginzburg-Landau free energy, numerically solved by the Simulated Annealing method. 



REFERENCES 

^M. Baert et al., Phys. Rev. Lett. 74, 3269 (1995). 

2 J.-Y. Lin et al., Phys. Rev. B 54, R12717 (1996). 

3 J.L Martin et al., Phys. Rev. Lett. 79, 1929 (1997). 

^G.S. Mkrtchyan, and V.V. Schmidt, Soviet Physics JETP 34, 195 (1972). 

5 A.I. Buzdin, Phys. Rev. B 47, 11416 (1992). 

^D. Saint- James and P.G. de Gennes, Phys. Lett. 7, 306 (1963). 

^S.M. Maurer, N.-C. Yeh, and T.A. Tombrello, Phys. Rev. B 54, 15372 (1996). 

^A. Bezryadin, A. Buzdin, and B. Pannetier, Phys. Lett. A 195, 373 (1994); Phys. Rev. B 
51, 3718 (1995). 

9G.M. Braverman, S.A. Gredeskul, and Y. Avishai, Phys. Rev. B 57, 13899 (1998). 

10 S. Kirkpatrick, CD. Gelatt, Jr., and M.P. Vecchi, Science 220, 671 (1983). 

i^M.M. Doria, J.E. Gubernatis, and D. Rainer, Phys. Rev. B 41, 6335 (1990). 

i^M.M. Doria, J.E. Gubernatis, and D. Rainer, Phys. Rev. B 39, 9573 (1989). 

i^Mauro M. Doria and Sarah C.B. de Andrade, Phys. Rev. B 53 3440 (1996). 

i^U. Klein and B. Pdttinger, Phys. Rev. B 44, 7704 (1991). 

""^^ In the discrete space integration becomes J{dv/V) -^ {'^/Np)J2n^=iJ2n^=i- Discretized 
gauge invariant derivatives are given in these reduced units by DyA -^ (A' — 
exp{2TTiAya/^o)A)/a, and D^A -^ {A" — A) /a. Our choice of gauge is {A^ = Q,Ay = 
(N^q/JL) Tix). In the kinetic energy the functions (r + r')/2 and (r + t")/2 are used, 
instead of just r, so to keep the symmetry with respect to nearest neighbor points. The 
functions (r',A') and (r",A") refer to the upper and right nearest neighbor sites of (t,A). 



^^ The present Simulated Annealing calculation was done with a typical value of 1, 600 Monte 
Carlo visits to each site at a given temperature reduction. 



FIGURES 
FIG. 1. The M vs. B curve is shown for the parameters R = 2,^0) J = 30, L = 20^0) T = 0, 

and K = 5. The number of FL in the unit cell varies from A^ = 1 to iV = 85. The fields Bp mark 

the transition to P trapped FL inside the CD. The inset shows the transition at iV = 68. 

FIG. 2. The Finh vs. B curve splits into five independent lines each one associated to P trapped 
FL inside the CD. The same parameters of Fig.[l] are used here. The inset shows the same quantity 
for the higher temperature of T/Tc = 0.6. 

FIG. 3. The two data sets correspond to the onset of localized superconductivity (white circle) 
and its disappearance (black circle). The inset shows the transition for B = 0.377, which occurs 
for T/Tc = 0.631. 

FIG. 4. The local supercharge {us = |Ap) and the local kinetic energy 

(ri)^ = [ 1 — r(x) ]^(T) |Z)Ap) densities are shown here for the same B {N = 24), but for two 
different temperature t = T/Tc, t = 0.615 and t = 0.630, below and very close the transition to the 
localized superconducting state. 
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